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SUMMARY 

The  development  of  a  turbulent  boundary  layer  flow  through  a  gap  in  an  isolated 
wall-mounted  roughness  element  has  been  studied  experimentally.  T wo  flow  regions  were 
distinguished  downstream  of  the  gap;  a  distortion  region  followed  b  a  readjustment  region. 
In  the  distortion  region  two  counteraeting  distortion  mechanisn.s  were  identified,  the 
relative  importance  of  which  depended  on  gap  size.  Thus  flows  downstream  of  large  gaps 
were  found  to  differ  significantly  from  those  through  small  gaps.  After  distortion  the 
layer  readjusts  itself  and  approaches  equilibrium  conditions  of  an  undisturbed  zero  pressure 
gradient  layer.  The  readjustment  starts  near  the  wall  with  the  turbulence  adjustment 
preceding  the  mean  flow  adjustment.  The  growth  of  the  internal  layer  for  flow  through  six 
different  gap  sizes  can  be  described  by  a  single  function  if  internal  layer  height  and  distance 
from  the  gap  are  non-dimensionalized  with  the  local  wall  length  scale.  Well  downstream  of 
the  gap  it  is  shown  that  all  six  flows  are  similar  and  arc  approaching  equilibrium  con¬ 
ditions  in  a  similar  manner. 
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NOTATION 


A  Constant  in  logarithmic  law  of  the  wall  (5.1) 

ci'  Local  skin  friction  coefficient  (to/£p«i2) 

cro'  Skin  friction  coefficient  of  undisturbed  layer 

D  Gap  width 

G  Clauser's  profile  shape  factor 

H  Height  of  roughness  element  (8  •  38  mm) 

u  Mean  velocity  in  X  direction 

^  u' a  Longitudinal  turbulence  intensity 

\/'um  2  Maximum  value  of  \/ u "i 
ur  Friction  velocity  (\  r0lp) 

ur(l  Friction  velocity  of  undisturbed  layer 

Ui  Free  stream  velocity  in  X  direction 

x  Distance  in  the  main  flow  direction 

i  Vertical  distance  from  wail 

ym'  Value  of  v  at  \  W2 

r  Local  wall  length  scale 

8  Boundary  layer  total  thickness 

8i  Height  of  internal  layer 

k  von  Karman  constant  (0-41) 

v  Kinematic  viscosity  of  fluid 

P  density  of  fluid 

m  Wall  shear  stress 


1.  INTRODUCTION 


There  has  been  considerable  work  on  the  response  of  turbulent  boundary  layers  to  various 
perturbations,  including  changes  in  wall  roughness  and  pressure  gradients,  changes  in  wall 
suction  and  injection  and  to  obstacles  in  the  flow.  Tani20  and  Schofield14  have  reviewed  this 
work  which  has  both  theoretical  and  practical  applications.  The  theoretical  interest  comes  from 
Clauser's7  proposal  that  these  experiments  offer  the  possibility  of  understanding  the  internal 
turbulent  exchange  mechanisms  in  the  boundary  layer  by  observing  the  response  of  a  layer  to 
a  step  or  pulse  input  in  one  of  its  boundary  conditions.  Practical  interest  in  the  experiments  arises 
from  the  need  to  predict  engineering  flows  involving  these  perturbations,  for  example  in  the  fields 
of  heat  transfer,  aerodynamics  of  jet  engines,  meteorology  and  wind  loads  on  buildings. 

The  present  work  considers  flow  of  a  turbulent  boundary  layer  over  and  through  a  gap  in 
a  surface-mounted  roughness  element,  which  has  practical  significance  for  example,  in  describing 
the  effects  of  buildings  and  other  obstructions  on  wind  gradients  and  turbulence  levels  near 
V/STOL  aircraft  landing  zones.  Data  on  such  flows  have  not  been  previously  reported,  but 
there  has  been  work  on  isolated  surface  roughness  elements  without  gaps.  Mueller  and 
Robertson11  studied  a  flat  plate  flow  over  a  wedge-shaped  element  the  height  of  which  was  of 
the  same  order  as  the  boundary  layer  thickness,  and  measured  mean  and  fluctuating  velocity 
profiles  up  to  40  element  thicknesses  downstream  of  the  element.  Reattachment  of  the  flow 
separation  caused  by  the  element,  defined  as  the  position  where  ct  —  0,  occurred  at  seven 
element  heights  downstream  of  the  element,  after  which  the  skin  friction  was  found  to  increase 
rapidly  to  its  undisturbed  value.  The  mean  and  fluctuating  velocity  profiles  were  judged  to  have 
returned  to  their  undisturbed  forms  at  a  distance  of  40  element  heights  downstream  of  the 
element.  Oka  and  Kostic12  performed  a  similar  study  using  a  square  roughness  element  on  the 
floor  of  a  rectangular  channel,  in  which  reattachment  was  again  found  to  be  located  at  seven 
element  heights  downstream  of  the  element.  Turbulence  intensity  profiles  after  reattachment 
showed  a  large  maximum  that  decreased  in  magnitude  and  moved  away  from  the  wall  with 
distance  downstream.  However,  these  maximum  values  at  18  element  heights  downstream  of 
the  element,  were  still  considerably  larger  than  values  in  an  undisturbed  layer.  Similar  results 
were  reported  by  Counihan,  Hunt  and  Jackson5  who  also  used  a  square  cross-section  element. 
These  authors  presented  an  eddy  viscosity  theory  that  adequately  described  the  mean  flow 
profiles  but  gave  poor  predictions  for  the  shear  stress  and  turbulence  intensity  profiles.  Results 
by  Sami  and  Liu15  suggest  that  the  distance  to  reattachment  behind  an  element  whose  height  is 
the  same  as  the  boundary  layer,  depends  on  the  Reynolds  number.  However,  their  reattach¬ 
ment  lengths  were  much  larger  (12  to  14  element  heights)  than  those  given  in  other  papers. 
Phataraphruk  and  Logan14  found  that  the  reattachment  length  behind  a  single  roughness 
element  in  a  smooth  pipe  depended  on  the  height  to  width  ratio  of  the  roughness  element,  which 
may  explain  the  different  reattachment  lengths  reported  by  authors  for  flows  behind  thin  two- 
dimensional  fences  (Petryk  and  Brundrett,13  17  fence  heights;  Good  and  Joubert,*  12-14  fence 
heights). 

Some  measurements  in  flows  with  no  gap  were  included  in  the  present  study.  Not  only 
does  this  two-dimensional  element  case  represent  a  limit  for  the  experiments  with  a  gap  but  by 
comparison  with  the  results  of  previous  work  these  tests  gave  confidence  in  the  results  for  flows 
through  gaps. 

Of  less  relevance  to  the  present  study  is  the  considerable  amount  of  work  which  has  been 
done  on  flow  around  isolated  three-dimensional  elements,  such  as  measurements  around  a 
surface-mounted  cube  by  Castro  and  Robins,4  around  model  buildings  by  Sharan,1*  around 
a  right  cylinder  by  Homung  and  Joubert*  and  around  a  surface-mounted  aerofoil  by  Brown.* 
These  flows  all  generate  horseshoe  vortices  around  the  front  of  the  models  that  shed  from 
both  sides  of  the  model  and  decay  fairly  rapidly  with  distance  downstream.  The  relevance  of 
these  studies  to  the  present  work  is  that  in  flow  through  a  gap,  vortices  are  shed  from  each  side 
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of  the  gap  and  have  a  significant  influence  on  flow  development  in  the  region  immediately  down¬ 
stream  of  the  gap. 


2.  EXPERIMENT 

A  two-dimensional  turbulent  boundary  layer  was  generated  and  its  (centreline)  develop¬ 
ment  observed  through  a  gap  in  a  wall-mounted  roughness  element  as  shown  schematically  in 
Figure  1.  The  element  was  square  in  cross-section  (8-38  mm  by  8  38  mm)  and  the  centreline 
gap  width  was  set  at  six  different  values  ranging  from  4- 19 mm  to  83-8 mm  (i.e.  0-5 H  to 
10  OH).  The  experimental  flow  was  generated  on  the  floor  of  a  suction  wind  tunnel  7  - 32  m  long 
and  0-56  m  wide.  The  flexible  roof  of  the  wind  tunnel  was  adjusted  to  give  a  zero  pressure 
gradient  flow  over  the  length  of  the  wind  tunnel.  Air  for  the  tunnel  was  firstly  filtered  and  then 
entered  a  honeycomb  flow  straightener.  i  50  mm  downstream  of  the  flow  straightener  the 
boundary  layer  was  tripped  by  a  12-7  mm  quaner  round  rod  and  then  further  thickened  by 
0-6  m  of  rough  sand  paper  fixed  to  the  floor.  A  smooth  wall  boundary  layer  was  allowed  to 
develop  over  a  length  of  4  m  upstream  of  the  roughness  element  under  study,  leaving  2-4  m 
of  tunnel  floor  for  observation  of  the  boundary  layer’s  response  to  the  element.  The  Reynolds 
number,  based  on  the  development  length  of  the  flow  at  the  element,  was  2-2x  10s,  which  was 
held  constant  for  all  tests.  The  total  boundary  layer  thickness  immediately  before  the  element 
was  77  mm  (~9 H). 

Pitot  tubes  and  hot  wire  anemometer  probes  were  mounted  on  a  traversing  mechanism 
above  the  roof  of  the  tunnel,  the  probes  entering  the  flow  through  slots  in  the  roof.  The  traversing 
mechanism  could  be  moved  only  along  the  tunnel  centreline  and  hence  only  profiles  along  the 
centreline  of  the  flow  were  recorded. 

As  a  test  of  the  flow  and  the  measuring  system,  mean  velocity  and  longitudinal  turbulence 
intensity  profiles  were  measured  in  an  empty  tunnel.  These  profiles  were  compared  with  zero 
pressure  gradient  profiles  by  Wieghardt  and  Tillmann22  and  Klebanoff  and  Diehl,10  and  gave 
satisfactory  agreement. 


3.  RESULTS 
3.1  Flow  Description 

As  only  centreline  profiles  were  measured  in  this  work  the  following  description  of  the  flow 
is  rather  speculative,  although  its  major  features  will  be  fairly  uncontroversial  at  least  in  quali¬ 
tative  terms. 

Three  regions  can  be  distinguished  in  the  flow  which  is  shown  schematically  in  Figure  I. 
The  first  region  is  well  upstream  of  the  element,  where  there  is  an  undisturbed  zero  pressure 
gradient  boundary  layer.  In  the  second  flow  region  the  layer  flows  through  the  gap  and  is  dis¬ 
torted  by  it.  The  distortion  starts  before  the  gap  and  continues  downstream  of  it  where  three- 
dimensional  separations  generated  at  each  gap  edge  bow  outwards  towards  the  flow  centreline. 
Thus  maximum  distortion  is  produced  at  a  position  a  little  downstream  of  the  element,  after 
which  the  distortion  rapidly  decreases  and  the  gap  edge  separations  join  up  with  the  separations 
caused  by  flow  over  the  roughness  element  (see  Fig.  1).  Results  suggest  that  most  of  the  dis¬ 
tortion  has  taken  place  at  the  longitudinal  position  where  the  separation  bubble  formed  by  flow 
over  the  roughness  element  closes.  Downstream  of  this  position  is  the  third  region  in  which  the 
layer  readjusts  itself  and  approaches  the  original  zero  pressure  gradient  form  which  is  its  final 
equilibrium  condition. 


3.2  Distortion  Region 

Distortion  of  the  centreline  flow  through  the  gap  is  caused  by  flow  adjacent  to  the  centreline 
funnelling  through  the  gap.  Here  the  effective  lateral  boundaries  to  the  flow  are  formed  by  sep¬ 
arated  regions  shedding  from  the  gap  edges,  almost  certainly  involving  counter-rotating  helical 
vortices  on  each  side  of  the  centreline  flow.  This  flow  picture  can  explain  the  different  behaviour 
of  the  centreline  flow  for  large  and  small  gap  widths.  Flow  through  the  large  gaps  is  well  removed 
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from  the  gap  edges  and  is  little  affected  by  the  gap  edge  separations.  Thus  the  flow  funnels 
through  the  gap,  increasing  the  mean  velocity  of  the  centreline  flow  and  decreasing  the  longi¬ 
tudinal  turbulence  intensity  as  in  a  conventional  nozzle.  On  the  other  hand  the  centreline  flow 
through  a  small  gap  is  affected  by  the  gap  edge  separations  which  in  this  case  occupy  a  significant 
proportion  of  the  flow  area  of  the  gap.  At  the  smallest  gap  settings  considered  here  the  gap 
edge  separation  vortices  dominate  the  gap  flow,  the  centreline  flow  is  brought  nearly  to  separation 
and  much  of  the  “through  gap”  flow  is  diverted  up  and  over  the  gap,  increasing  mean  flow 
velocities  at  distances  greater  than  the  element  height.  The  gap  edge  separations  adjacent  to  the 
centreline  flow  increase  its  turbulence  level  near  the  wall.  Data  supporting  this  flow  picture 
are  presented  in  Figure  2  where  a  series  of  profiles  for  a  small  gap  ( D/H  =10)  and  for  a  large 
gap  (D/H  =  5  0)  are  compared  with  the  (upstream)  undisturbed  boundary  layer  profiles.  In 
the  distortion  regime  (0  <  x/H  <  10)  the  small  gap  profiles  show  a  reduction  in  mean  velocity 
near  the  wall  and  the  corresponding  profiles  of  longitudinal  turbulence  (\/u'*juro)  show  a  large 
increase  near  the  wall.  Note  that  both  these  series  of  profiles  show  maximum  distortions  at 
about  four  to  six  element  heights  downstream  of  the  gap.  This  is  due  to  the  bowing  out  of  the 
separation  boundaries  downstream  of  the  gap. 

The  profiles  for  the  large  gap  show  the  opposite  behaviour.  In  the  distortion  regime  the 
gap  edge  separations  are  well  removed  from  the  centreline  and  thus  mean  flow  velocities  are 
increased  (relative  to  the  undisturbed  profile)  and  longitudinal  turbulence  intensities  decreased. 
Note  that  in  this  case  maximum  distortion  occurs  in  profiles  nearest  the  gapt  as  it  is  the  actual 
gap  not  the  downstream  separation  boundaries  that  is  causing  the  majority  of  the  centreline 
distortion. 

In  both  cases  the  outer  portions  ( i  ,  H  >  2)  of  the  profiles  at  the  end  of  distortion 
(x/H  ~  10)  are  again  coincident  with  the  undisturbed  upstream  profiles.  This  is  fortunate  as 
it  aids  the  analysis  of  the  layer's  readjustment  that  is  discussed  in  the  next  section. 

The  existence  of  vortices  shed  from  the  gap  is  supported  by  centreline  profiles  taken  in  the 
smallest  gap  considered  ( D/H  0 - 5).  The  wall  regions  of  mean  and  fluctuating  velocity  profiles 
at  x/H  =  0  are  --hown  in  Figure  3.  Corresponding  positions  in  the  two  profiles  (i.e.  the  same 
values  of  y/ H )  arc  marked  by  circled  numbers  at  salient  profile  points.  Figure  3 a  shows  the  mean 
profile  plotted  on  Clauser's**  semitogarithmic  coordinates,  u,U\  versus  \ogu>  y(Ui/v).  For  a 
normal  two-dimensional  boundary  layer  this  mean  flow  profile  should  show  a  linear  distribution 
up  to  about  0158.  However,  in  this  case  the  profile  shows  pronounced  sinuosity.  The  peaks 
and  valleys  of  this  sinuosity  are  closely  related  to  sinuosity  in  the  longitudinal  turbulence 
intensity  profile  (Fig.  3b).  The  sinuosities  can  be  followed  in  downstream  profiles  (for  this  gap 
setting)  but  show  a  rapid  decay.  These  data  give  good  support  to  the  postulate  of  significant 
gap  edge  separations  and  their  associated  vortices  affecting  the  centreline  flow  at  small  gap 
settings. 


3.3  Readjustment  Region 

The  readjustment  of  the  centreline  flow  towards  velocity  distributions  of  an  undisturbed 
zero  pressure  gradient  layer  start  at  about  x/H  10  for  all  gap  settings.  At  this  position  dis¬ 
tortion  of  the  profiles  by  the  gap  is  small  and  rapidly  decaying.*  Downstream  of  x/H  ~  10  the 
dominant  feature  of  the  flow  is  readjustment,  which  under  these  conditions  is  in  many  ways 
similar  to  the  adjustment  of  a  boundary  layer  to  a  new  wall  roughness  after  a  rough  to  smooth 
step  change  in  surface  roughness  as  studied  by  Antomu  and  Luxton.1  The  “new"  flow  con¬ 
dition  in  the  present  case  is  the  cessation  of  flow  distorting  forces  rather  than  the  cessation  of  a 
rough  wall.  It  is  well  established  (Tani;20  Schofield  Townsend'-')  that  flow  modification  occurs 
outwards  from  the  wall  in  all  readjusting  layers  as  the  scale  of  turbulence  near  the  wall  is  small, 
short  lived  and  hence  can  make  thi  most  rapid  adjustment  to  the  new  flow  conditions.  Further 
out  from  the  wall,  as  well  as  small  eddies  there  exist  large  eddies,  which  contain  most  of  the 


t  In  Figure  2  this  is  the  profile  a!  x/H  2.  Earlier  profiles  at  x/H  ---  I.  not  shown  in 
Figure  2,  have  slightly  more  distortion. 

+  This  can  be  seen  to  a  limited  extent  in  Figure  2  but  better  evidence  is  presented  later  in 
connection  with  Figure  7. 


turbulent  energy,  have  longer  lifetimes  and  travel  downstream  faster  than  wall  eddies  because 
they  exist  in  regions  of  higher  mean  velocity  (Townsend21).  For  these  reasons  the  readjustment 
process  proceeds  ever  more  slowly  as  flow  modification  works  its  way  outwards  through  the 
layer,  and  complete  attainment  of  the  final  equilibrium  profile  appears  to  take  a  very  long 
development  distance. t 

To  find  the  height  of  modified  flow  (usually  termed  the  height  of  the  internal  layer)  Si,  the 
readjusting  profile  is  compared  with  the  profile  the  layer  would  have  had  if  the  perturbation  had 
not  been  applied  to  the  layer  (see  Schofield18,17).  In  the  present  case  the  unperturbed  profile  is 
simply  that  of  a  zero  pressure  gradient  layer  and  this  is  compared  with  the  actual  profiles  in 
Figure  2.  The  rates  of  adjustment  of  various  flow  parameters  (such  as  mean  velocity,  fluctuating 
velocities,  turbulent  shear  stress)  are  not  necessarily  the  same.  However,  as  all  workers  in  this 
field  have  measured  mean  velocities  it  is  usual  to  determine  the  8|  associated  with  mean  velocity 
adjustment.  Thus  the  height  of  the  internal  layer  8!  is  defined  as  the  distance  from  the  wall  that 
the  actual  mean  velocity  profile  merges  with  the  unperturbed  profile,  as  shown  in  Figure  2.  It 
can  be  seen  that  modification  of  the  longitudinal  turbulence  intensity  profiles  precedes  modific¬ 
ation  of  the  mean  velocity  for  both  small  and  large  gap  cases.  This  result  applied  to  all  six  cases 
studied  in  this  work. 

Many  authors  have  expressed  the  variation  of  internal  layer  height  with  distance  from  the 
perturbation  in  a  straight  dimensional  form;  e  g.  Si  oc  43  by  Antonia  and  Luxton1  for  a  rough 
to  smooth  wall  step.  Such  dimensional  relations  may  not  be  useful  for  other  flows  with  different 
boundary  conditions,  geometries  or  Reynolds  numbers.  Schofield17,18  has  shown  that  non- 
dimensionalizing  the  results  with  the  local  wall  length  scale  (z)  introduced  by  Townsend,21 
collapses  data  from  a  very  wide  variety  of  flows  over  a  step  change  in  roughness  onto  a  single 
curve.  The  wall  length  scale  for  smooth  wall  flow  is  related  to  the  wall  shear  and  its  definition 
is  derived  from  the  logarithmic  law  of  the  wall  which  applies  to  the  inner  regions  of  most  two- 
dimensional  turbulent  boundary  layers.  It  can  be  written 


u/w,  =-•  (]/«•)  log tyiujv)  +  A 

(1) 

and  the  local  wall  length  scale  (z)  is  defined  by  rewriting  it  as 

U/UT  =  (I /k)  log«  y/Z 

(2) 

which  by  comparison  gives 

z  =  v/(uT  exp  M))  { 

(3) 

Using  the  local  value  of  z  for  each  profile  the  variation  of  the  non-dimensional  internal  layer 
height  (S|/z)  is  compared  with  the  non-dimensional  distance  from  the  gap  (x/z)  in  Figure  4.  It  is 
seen  that  the  data  expressed  in  this  form  are  similar,  although  a  systematic  variation  with  gap 
size  can  be  discerned.  This  systematic  variation  could  be  attributable  to  a  continuing  influence 
of  the  pulse  perturbation  that  varies  with  gap  size.  However,  all  the  data  collapse  onto  a  single 
curve  with  little  scatter  as  shown  by  the  last  graph  in  Figure  4. 

Townsend’s  analysis  considered  thick  zero  pressure  gradient  flow  over  a  small  step  change 
n  roughness  and  gave  the  following  relation  for  the  internal  layer  height 

(8,/z)(loge(8l/z)  -  1 )  =  2*2  (x/z).  (4) 

Scholfied17, 18  showed  that  this  equation  gave  quite  good  descriptions  for  flow  over  a  wide  range 
of  step  changes  in  roughness  under  all  types  of  flow  conditions,  many  of  which  were  very  different 
from  those  of  the  Townsend  analysis.  Equation  (4)  has  been  plotted  in  Figure  4  and  appears  to 
give  a  useful  description  of  the  internal  layer  growth  for  all  six  flows  even  though  the  present 
perturbation  was  of  a  different  kind  from  that  analysed  by  Townsend. 

Figure  5  compares  mean  velocity  and  longitudinal  turbulence  intensity  profiles  at  the  last 
measuring  station  ( x/H  =  224)  for  all  gap  sizes.  The  mean  profiles  for  all  six  cases  are  closely 


tThis  is  a  common  conclusion  of  Bradshaw  and  Wong,2  Townsend21  and  Schofield18,17,18 
and  is  further  demonstrated  later  in  this  paper. 

t  Values  of  uT  were  derived  by  the  Clauser8  method  which  is  discussed  later  in  this  paper. 
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similar  while  the  turbulence  intensity  profiles  are  less  so.  The  two  sets  of  profiles  for  the  previous 
station  ( x/H  =112)  are  also  similar  but  the  degree  of  similarity  is  less  marked.  Thus  it  appears 
that  during  the  readjustment  process  these  six  layers  have  “forgotten”  the  difference  in  their 
distortion  histories  and  tend  to  common  velocity  distributions  as  they  approach  their  final  equi¬ 
librium  form.  Figure  5  also  shows  that  both  mean  and  fluctuating  flows  have,  by  x/H  =  224, 
returned  to  their  equilibrium  distributions  in  the  region  immediately  adjacent  to  the  wall. 
However,  both  profiles  still  have  considerable  adjustment  to  make  before  they  are  identical  to 
the  undisturbed  zero  pressure  gradient  profiles.  Thus  these  results  do  not  agree  with  the  con¬ 
clusion  of  Mueller  and  Robertson11  that  the  flow  (for  the  zero  gap  case)  returns  to  undisturbed 
conditions  at  40  element  heights  downstream  of  the  element.  The  present  results  are  consistent 
with  Bradshaw  and  Wong’s2  conclusion  that  full  readjustment  takes  hundreds  of  boundary  layer 
lengths  downstream  of  an  obstacle. 


3.4  Mean  Flow  Parameters 


The  skin  friction  coefficients  (ct')  for  all  profiles  were  determined  by  plotting  the  mean 
velocity  data  on  axes  u/Ui  versus  logioy((7i/v)  and  comparing  the  wall  region  of  the  profile 
with  the  law  of  the  wall  (Equation  (1))  expressed  in  Clauser’s6  form: 

(u/U ,)  =  5 •76(ef72)1'2  logio )tUj/y)  +  5-76(r,'/2)1  2  log10  (r,72)1'2  +  5-\(c,'/2)  (5) 

Most  profiles  showed  good  agreement  with  the  Clauser  lines  and  hence  reliable  values  of  skin 
friction  coefficient  could  be  inferred.  Early  profiles  in  the  small  gap  flows  did  show  sinuousity 
as  was  discussed  with  reference  to  Figure  3  and  thus  the  skin  friction  coefficients  for  these  profiles 
were  more  doubtful. 

Ratios  of  local  skin  friction  coefficient  to  the  skin  friction  coefficients  of  an  undisturbed 
layer  are  plotted  in  Figure  6  for  all  gap  sizes  as  a  function  of  distance  from  the  gap.  During  the 
distortion  phase  (0  <  x/H  <  10)  different  wall  shear  distributions  are  produced  for  flows 
through  different  gap  sizes.  For  flow  through  small  gaps  ( D/H  <  5)  the  centreline  flow  approaches 
separation  under  the  influence  of  the  gap  edge  vortices,  so  that  the  skin  friction  coefficients 
approach  zero.  For  the  smallest  gap  studied  (D/H  —  0 - 5)  the  flow  only  just  remains  attached. 
As  gap  size  increases,  the  gap  edge  separations  are  moved  away  from  the  flow  centreline  and 
occupy  a  smaller  proportion  of  the  flow  area  in  the  gap;  thus  the  minimum  skin  friction  on  the 
centreline  increases.  At  the  largest  gap  size  (D/H  =  10  0)  the  increased  velocity  through  the  gap 
results  in  a  centreline  skin  friction  coefficient  that  exceeds  the  undisturbed  value  for  x/H  <  20. 

During  the  initial  readjustment  phase  the  skin  friction  distributions  depend  on  the  value 
at  the  end  of  the  distortion  phase.  However,  all  distributions  are  approaching  the  final  equi¬ 
librium  value  of  skin  friction  from  below  at  the  last  measuring  station.  The  behaviour  of  the 
skin  friction  distributions  at  x/H  =  224  is  not  however  uniform:  extrapolated  distributions  for 
small  gap  settings  overshoot  the  equilibrium  value  a  short  distance  downstream  of  x/H  <  224, 
whereas  those  for  the  two  largest  settings  do  not. 

For  the  zero  gap  case  Bradshaw  and  Wong2  also  reported  overshooting  of  the  equilibrium 
skin  friction  coefficient.  Another  point  of  agreement  of  these  results  with  previous  work  is  the 
position  of  reattachment  for  the  zero  gap  case  which,  as  shown  in  Figure  6,  occurs  at  approxi¬ 
mately  x/H  =  8.  This  broadly  agrees  with  previous  results  by  Mueller  and  Robertson11  (reattach¬ 
ment  at  x/H  ~  7),  Oka  and  Kostic12  (at  x/H  ~  7)  and  Counihan  et  a/.s  (at  x/H  ~  6). 

An  accurate  mean  velocity  profile  shape  factor  was  introduced  by  Clauser®  and  is  defined  as 


c=irr 


(6) 


For  an  undisturbed  zero  pressure  gradient  layer  at  high  Reynolds  number,  G  has  a  value  of  about 
6-8.  Perturbations  applied  to  a  zero  pressure  gradient  layer  cause  G  to  vary  from  this  equi¬ 
librium  value.  The  variation  of  G  after  a  perturbation  is  a  useful  indication  of  the  degree  of 
distortion  and  subsequent  readjustment  of  the  whole  layer.  Distributions  of  G  for  all  gap  sizes 
in  the  present  experiment  have  been  plotted  in  Figure  7.  Not  surprisingly,  maximum  distortion 
(i.e.  largest  G  variation)  occurs  downstream  of  the  smallest  gaps.  However,  the  profile  distortion 
is  remarkably  similar  for  all  gap  sizes  after  about  x/H  -■  60,  which  reinforces  the  previous 
conclusion  that  layers  forget  their  distortion  history  in  the  readjustment  region  and  follow 
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similar  paths  back  to  equilibrium  conditions.  There  are,  however,  minor  differences  in  G  vari¬ 
ation  at  large  distances  from  the  gap.  Values  of  G  for  small  gap  flows  undershoot  the  equilibrium 
value  at  xjH  =  224  whereas  for  the  largest  gaps  G  approaches  the  equilibrium  value 
monotonically.  Undershooting  for  the  zero  gap  case  has  been  previously  reported  by 
Bradshaw  and  Wong,2  who  found  that  the  minimum  value  of  G  occurred  at 

xjH  =  I00VS/W 

which  for  the  present  experiment  gives  x/H  —  316  as  the  position  for  minimum  G  which  seems 
reasonable  for  the  data  presented  in  Figure  7. 


3.5  Turbulent  Flow 

Figure  8  plots  the  maximum  non-dimensional  longitudinal  turbulent  intensity  (v  um'2/uro) 
for  each  profile,  and  the  height  of  this  maximum  from  the  wall,  as  a  function  of  distance  from  the 
gap.  There  are  again  different  behaviours  for  flow  through  small  and  large  gaps.  Small  gaps 
produce  large  maximum  values  of  \  «m'~  which  are  attained  shortly  after  the  gap  and  are 
(approximately)  maintained  throughout  the  distortion  region,  vum'2  decays  from  this  maximum 
value  in  the  readjustment  region.  The  graphs  are  truncated  at  x/H  =  56  because  flow  readjustment 
produces  profiles,  downstream  of  x/H  =  56,  in  which  y  um‘-  reverts  to  a  position  very  near  the 
wall  as  in  the  undisturbed  upstream  profile  (sec  Fig.  2). 

Figure  8  also  plots  vm'  (the  distance  from  the  wall  to  \  um'2)  which  is  initially  constant 
and  then  suddenly  increases  with  increasing  distance  downstream  of  the  gap.  This  sudden  increase 
commences  at  progressively  larger  distances  as  the  gap  is  increased.  For  the  largest  gap 
(D/H  —  10  0)  an  increase  in  ym'  does  not  occur  within  the  observation  range. 

These  results  are  consistent  with  flow  mechanisms  previously  described.  Increasing  values 
of  rm'  with  distance  reflect  a  flow  adjustment  that  works  outwards  from  the  wall.  High  initial 
values  of  \  Um'-  for  small  gap  flows  reflect  the  higher  turbulence  levels  generated  by  the  gap  edge 
separations.  For  large  gap  flows,  turbulence  is  generated  at  the  beginning  of.  and  as  part  of,  the 
readjustment  process.  These  large  gap  turbulence  profiles  thus  have  a  small  maximum  which 
appears  near  the  end  of  distortion  or  at  the  beginning  of  readjustment. 


4.  CONCLUSIONS 

(1)  Boundary  layer  flow  through  a  gap  in  a  small  isolated  roughness  element  consists  of  a 
flow  distortion  region  followed  by  a  flow  readjustment  region.  The  distortion  region 
extends  about  10-16  element  heights  downstream  of  the  element.  The  following 
readjustment  region  appears  to  extend  for  several  hundred  element  heights  downstream 
of  the  element. 

(2)  Centreline  flow  through  small  gaps  has,  compared  with  an  undisturbed  flow: 

(i)  an  increased  longitudinal  turbulence  intensity  near  the  wall; 

(ii)  a  decreased  mean  velocity  near  the  wall; 

(iii)  a  large  decrease  in  wall  shear  up  to  x/H  —  50:  and 

(iv)  a  large  change  in  mean  profile  shape  up  to  xjH  -  50. 

Conversely,  the  centreline  flow  through  large  gaps  has 

(i)  a  decreased  longitudinal  turbulence  intensity  near  the  wall; 

(ii)  an  increased  mean  velocity  in  the  distortion  region;  and 

(iii)  little  change  in  wall  shear  or  in  mean  profile  shape. 

(3)  The  mechanism  of  flow  readjustment  downstream  of  distortion  proceeds  by  a  flow 
modification  that  works  outwards  from  the  wall.  Both  the  fluctuating  and  mean  flow 
fields  arc  modified;  the  modification  of  the  fluctuating  flow  precedes  the  mean  flow 
modification  for  ail  gap  sizes.  The  growth  of  this  new  internal  layer  as  a  function  of 
distance  from  the  gap  can  be  collapsed  (for  all  six  cases)  onto  a  single  curve  with  little 
scatter  if  the  parameters  are  ncn-dimensionalized  with  the  local  wall  length  scale.  This 
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curve  is  usefully  approximated  bv  a  relation  derived  by  Townsend  lor  How  over  a 
roughness  perturbation. 

(4)  Initial  readjustment  of  centreline  flow  parameters  (cf  and  (>)  depends  on  their  values 
at  the  end  of  distortion.  However,  iter  readjustment  behaviour  is  similar  for  all  gap 
sires  with  the  parameter  curves  and  the  mean  and  fluctuating  profiles  being  similar  for 
v  H  >  100. 

(5)  Boundary  layer  flow  over  an  isolated  element  without  a  gap  can  be  regarded  as  a  limiting 
case  of  boundary  layer  flow  through  a  gap.  In  the  zero  gap  case,  profile  shapes  and  para¬ 
meter  variations  downstream  of  reattachment  follow  the  pattern  set  by  centreline  flows 
through  small  gaps. 
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FIG.  1  SCHEMATIC  VIEW  OF  EXPERIMENTAL  FLOW 


SELECTED  MEAN  VELOCITY  AND  LONGITUDINAL  TURBULENCE  INTENSITY 


FIG.  3  MEAN  VELOCITY  AND  LONGITUDINAL  TURBULENCE  INTENSITY  PROFILES 
FOR  D/H  =  0.5  at  x/H  =  0 


VELOCITY  AND  LONGITUDINAL  TURBULENCE  PROFILES  AT  x/H  «  224 
LL  GAP  SIZES 
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FIG.  8  MAXIMUM  VALUES  OF  LONGITUDINAL  TURBULENCE  INTENSITY  AND 
DISTANCE  FROM  WALL  OF  THESE  MAXIMA 
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